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We consider the properties of a self-avoiding polymer chain, adsorbed on a solid attractive sub- 
strate which is attached with one end to a pulling force. The conformational properties of such 
chain and its phase behavior are treated within a Grand Canonical Ensemble (GCE) approach. 
We derive theoretical expressions for the mean size of loops, trains, and tails of an adsorbed chain 
under pulling as well as values for the universal exponents which describe their probability distri- 
^ ' bution functions. A central result of the theoretical analysis is the derivation of an expression for 

| the crossover exponent cj>, characterizing polymer adsorption at criticality, <j> = a — 1, which relates 

the precise value of 4> to the exponent a, describing polymer loop statistics. We demonstrate that 
1 — Til < a < 1 + u, depending on the possibility of a single loop to interact with neighboring loops 
in the adsorbed polymer. The universal surface loop exponent 711 « —0.39 and the Flory exponent 
v ?a 0.59. 

We present the adsorption-desorption phase diagram of a polymer chain under pulling and demon- 
strate that the relevant phase transformation becomes first order whereas in the absence of external 
force it is known to be a continuous one. The nature of this transformation turns to be dichotomic, 
i.e., coexistence of different phase states is not possible. These novel theoretical predictions are 
verified by means of extensive Monte Carlo simulations. 
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I. INTRODUCTION 



C3 

S 

With the development of novel single macromolecule experiments, the manipulation of individual polymer chains 
and biological macromolecules is becoming an important method for understanding their mechanical properties and 
characterizing the intermolecular interactions P, 0]. Much of the related upsurge of interest into the statics and 
dynamics of single macromolecules at surfaces has been spurred by the use of Atomic Force Microscopy 0, 0, H, @ 
(AFM) and optical/magnetic tweezers 0, d, HI which allow one to manipulate single polymer chains. Measurements 
of the force, needed to detach a chain from an adsorbing surface, and most notably, of the force versus extension 
relationship which exhibits sharp discontinuities have been interpreted as indication for the presence of unadsorbed 
loops on the surface. In turn, this has initiated a number of theoretical studies [13, HH, 03 which have helped to get 
better insight into the thermodynamic behavior and the mechanism of polymer detachment from adhesive surface 
under pulling external force. A comprehensive treatment of the problem for the case of a phantom polymer chain 
can be found in the paper of Skvortsov et al. [l3j]. There is a close analogy between the forced detachment of 
an adsorbed polymer chain like polyvinilamine and polyacrylic acid, adhering to a solid surface such as mica or a 
self- assembled monolayer, when the chain is pulled by the end monomer, and the unzipping of homogeneous double- 
stranded DNA. In the context of DNA denaturation and the simple single chain adsorption this analogy has been 
discussed already in the middle 60s fl4j . Recently, the DNA denaturation and its unzipping have been reconsidered by 
Kafri, Mukamel and Peliti [l5j]. The consideration was based on the Poland and Sheraga' s Grand Canonical Ensemble 
(GCE) approach Tfl 17| as well as on Duplantier's analysis of the number of configurations in polymer networks of 



arbitrary topology [la ]. Duplantier's analysis makes it possible to calculate the values of universal exponents which 
undergo renormalization due to excluded volume effects. In particular, it has been shown by Kafri et al. [H[ that this 
renormalization procedure changes even the order of the melting (or denaturation) transition in DNA from second to 
first order. 

In the present paper we use the approach of Kafri et al. [H| in order to treat the detachment of a single chain from 
a sticky substrate when the chain end is pulled by external force. It has been pointed out earlier [13j that the problem 
may be considered within the framework of two different statistical ensembles, i.e., by keeping the pulling force fixed 
while measuring the (fluctuating) position of the polymer chain end, or, by measuring the (fluctuating) force necessary 
to keep the chain end at fixed distance above the adsorbing plane. Our theoretical consideration has been carried out 
in the fixed force ensemble whereas experimentalists usually work in the fixed distance ensemble. We start in Section 
HI1 with the consideration of the conventional adsorption (i.e. force-free) problem where we derive a basic expression 
for the crossover exponent describing polymer adsorption. There we also consider theoretically some basic features of 
adsorbed polymer chains as the variation of the average length of loops and tails in the chain with changing strength 
of the adsorption potential. In Section [Ml we extend our theoretical analysis to the case of polymer adsorption in the 
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presence of external force, and obtain results for the main conformal properties of such chains as well as the relevant 
phase diagram of the system. The properties of the simulation model are briefly reviewed in Section HVi and then in 
Section [V] we report on our most important results, gained in the course of the computer experiment, and compare 
them to theoretical predictions. We end this work in Section I VII with a brief summary and discussion of the most 
salient results of the present investigation. 



II. SINGLE CHAIN ADSORPTION: LOOP-, TRAIN-, AND TAIL STATISTICS 

A single chain, adsorbed on a solid plane, is built up from loops, trains and a free tail. In order to derive expressions 
for the mean values of these basic structural units, one may treat the problem within the Grand Canonical Ensemble 
(GCE)[l7l. In the GCE approach the lengths of these building blocks are not fixed and are allowed to fluctuate. The 
GC-partition function is given as 



E(z) 



Vo(z) Q(z) 
1 - V(z)U(z) 



(2.1) 



where Sat is the canonical partition function of a chain of length N and z is the fugacity. U (z), V(z) and Q(z) denote 
the GC partition functions of loops, trains and a tail respectively. The building block adjacent to the tethered chain 
end is allowed for by Vq{z) = 1 + V(z). The series given by Eq. (|2.1|) is a geometric progression with respect to 
U(z)V(z). Figure [1] gives a pictorial representation of this series. 



Vn 



+ 





FIG. 1: Schematic representation of the series expansion given by Eq. (|2.1|l 



The GC-partition function of the loops is defined by 

oo 

U(z) = 



m l-7n 



m—1 



(2.2) 



where fi TO is the number of surface m-loops (i.e., self-avoiding walks of length m which start and terminate on the 
surface) configurations. For an isolated m-loop this number of configurations is given by fi m = [j,™m' ril ~ 1 where Ma 
is the 3c? connective constant (in three dimensions, d = 3, one has = 4.68, and the exponent 711 = —0.390 )[l9(. 
Below we will demonstrate that the exponent 711 changes due to the excluded volume interactions between different 
loops. 

The train GC-partition function reads 



V(z) = 



OO 

E 



(p,2Wz) m 

m. 1 1d = 2 



(2.3) 



where the number of train configurations of length m (which are located in the d = 2 surface plane) is given by 
* m = w m m 7 ^ 2-1 . Here fi 2 = 2.6 and ~f d=2 = 1.343 [H[. In Eq. (|2~3")1 we have taken into account that 
each adsorbed segment of the chain gains an additional statistical weight w = exp(e/fcsT) = exp(e), where T is 
the temperature and the Boltzmann constant ks is set to unity. In what follows the notation e stands for the 
dimensionless adsorption energy of a single monomer. In fact, e denotes the potential well depth of the short-ranged 
surface potential, defined in the description of our simulation model in Section IIVI 
The GC-partition function for the chain tail is given by 



q(*)=i+e A m z™=i+x: s 



m — 1 m—1 

where the m-tail number of configuration equals A m = fi™ m 7l_1 , and in d = 3 the exponent 71 = 0.680 [l| 



(2.4) 
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With the knowledge of the GC partition function, given by Eq. (|2.ip . it is possible to calculate the number of 
weighted configurations of a polymer chain, containing TV segments (i.e., its canonical partition function), Ejy. From 
the generating function method (see, e.g., Sec. 2.4 in the book by Rudnick and Gaspari [2(|) it is well known that at 
N — > oo the coefficient at z N is defined by a singular point (a pole or a branching point) of S(z) which lies closest to 



the origin. In our case this is a simple pole, 



which is determined from the condition 

V(z*) U{z*) = 1 



The principal contribution to this coefficient at z N is (z*) 
energy 



i.e., ^ N 



CO 



F 



-TlnS 



TNkiZ* 



(2.5) 

and so the corresponding free 

(2.6) 



In Section |TT1 devoted to the adsorption of a pulled polymer chain, we shall see that an important singularity arises 
also from the tail generating function. The average fraction of adsorbed monomers, n = N s /N (where N s is the 
number of adsorbed monomers) which we use as an order parameter for the degree of adsorption, can be calculated 
then as follows 



N 



1 d lnSjy 
N dlnw 



dhi, 



dhiw 



(2.7) 



The generating functions, given by Eqs. (|2.2p . (12. 3p . and (|2.4p , can be conveniently expressed in terms of the polylog 
function 21| . In the Appendix we sketch the properties of the polylog function and its behavior in the vicinity of the 
singular point. In terms of the polylog function (see Appendix) the basic Eq. (|2.5p is then given by 



$(a,/i 3 z*) = $ 1 (X,f Jl2 wz*) 



(2.8) 



where the exponents a = 1 — 711 rj 1.39 > 1 and A = 1 — jd=2 ~ —0.343 < 1. One should note that the exponent 
a = 1 — 711 corresponds to a loop treated as an isolated one. This is an important feature of the method which 
handles the main building blocks (loops, trains and tails) as independent objects (see, e.g., Eq. (|2.ip ). Nevertheless, 
in Sec. Ill B[ following Kafri et al. [la], we shall show that by taking into account the excluded volume interaction 
between a loop and the rest of the chain one ends up with a renormalized value of the exponent a (it increases). This 
is important because the value of a determines itself the value of the well known surface (or, crossover) exponent <j) 
in all the basic scaling laws pertaining to polymer adsorption (see below). 

Close to the critical point, z c = z* which is defined by ^z c = 1, the l.h.s. of Eq. (|2.8p can be expanded (cf. 
Ea. (|Alip ) as follows 

C(a) - a a (l - fi 3 z*) a ^ ~ Ml - M z l = * _1 (1 - 7d=2,»2wz*) (2.9) 

with C( x ) denoting the Riemann zeta-function. At the critical adsorption point (CAP), e c and w c = exp(e c ), the 
solution of Eq. (|2.8p is z* = z c — I//Z3 so that w c is given by the expression 

C{a) = <p- 1 (l- ld=2 ,fi 2 w c /v 3 )- (2.10) 

The expansion of Eq. (|2.9p around the critical point, z* — z c and w = w c , could be effected by the substitution of 
w = w c + 5 and z* = z c — A in Eq. (|2.9p . Here S and A are corresponding infinitesimal increments and we took into 
account that z* decreases with increasing w. Substituting this in Eq. (|2.9p gives 



C(a) -M/! 3 A) Q 1 0--ld=2,^2W c z c ) 



$-2(1 



ld=2,^2W c Z c ) 



dx 



6 (2.11) 



Taking into account the condition for the critical point, Eq. (|2.10p . as well as the identity Eq. 
z* can be recast in the form 



2l, the solution for 



z*(w) 



1 



l/(o-l) 



(2.12) 



where the constants 



.4 



M2$(-7rf=2, ^2W c /fl 3 ) 

$ 2 (1 - 7d=2,M2Wc/^3) 
7T 



r(a)| sin(7TQ!)| 



(2-13) 
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FIG. 2: The 'order parameter' (i.e., the fraction of adsorbed chain segments), n, against the surface potential, e, in the absence 
of detachment force, / = 0. The inset shows the variation of the fugacity z* with w — exp(e), Eq. (|2.8|l . 



and w c is denned by Eq. (|2.10p . The full numerical solution for the order parameter as well as for the pole z*(w) is 
displayed in Fig. [21 

Having the solution Eq. (|2.12| ) at hand, one can use the expression Eq. (|2.7|) for the average fraction of adsorbed 
monomers. After some straightforward calculations one arrives at 

n(e) oc (e-ej^r- 1 (2.14) 

where one has used w — w c k, exp(e c ) (e — e c ). On the other hand, it is well known [19j that the scaling behavior in the 
vicinity of the critical adsorption energy is described by the crossover exponent cf>. Namely, the corresponding scaling 
relationship is given by 

(n^- 1 , at e = e c , . 

n(e)oc \ L (2.15) 

[ (e - e c ) * at e > e c 

If the result, given by Eq. (|2.14p . is compared to that of Eq. (|2.15p . it becomes apparent that 

<j) = a-l (2.16) 

This result, derived first by Birshtein 17], is of principal importance. Here it is derived in the context of self-avoiding 
chains. As stated above, if the loops are treated as independent non-interacting objects, the exponent a = 1 — 711, 
so that 

= - 7ll w 0.39 (2.17) 

In Sec. IIIBI wc shall demonstrate that by taking into account the excluded volume interactions between a loop and 
the rest of the chain one finds an increase of the values of a, and <f>, respectively. 



A. Loops and tails distributions close to criticality 



Here we examine how the size distribution of polymer loops and tails looks like close to the critical point of 
adsorption. The GC-partition function for loops, given by Eq. (|2.2p . yields immediately 



Ploop(0 



(>3Z) 



*\l 



li+4 



(2.18) 



where we have used the essential relation between the loop exponent a and the crossover exponent cf>, Eq. (|2.16[ ). Close 
to the critical point, /13Z* < 1 (see Eq. (|2.12[ and the Z-dependence is mainly described by inverse power-law 



li+4> 



(2.19) 
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The power-law decay for the loop distribution close to the criticality has been discussed in the early 80s by P.-G.de 
Gennes [22| . Deeper in the region of adsorption, however, the exponential part in Eq. (|2.18[) dominates. Taking into 
account Eq. (|2.12[) . one obtains 

1 



-Ploop(0 



■ exp 



(2.20) 



i.e., with increasing adsorption energy e the size distribution becomes narrower. 

It is of interest to note that the size distribution of loops can be reformulated in terms of the distribution F(h) 
of projected lengths h of the loops between two consecutive monomers residing on the adsorbing surface which 
was analyzed by Bouchaud and Daoud [23}. The relation between Pi OO p(0 and F(h) is straightforward, namely 
F(h)dh — P\ oop (l)\dl / dh\dh, where due to the isotropy of loops h ~ V too. Taking these relations into account as 
well as Eq. (|2.19j) one obtains 



F(h) 



(2.21) 



This corresponds exactly to the result of Bouchaud and Daoud [23] where such broad distribution was associated with 
the so-called node- avoiding Levy flight. 

The distribution of tails (at the CAP, i.e., at fj, 3 z* < 1) is even broader, namely 



(2.22) 



where for an isolated tail /3 = 1 — 71 » 0.32. We will show below (see Eq. (l2.46p ) that if the interaction of a tail with 
the rest of the chain is taken into account this leads to a larger value of f3 = 0.51. One should be aware, however, 
that this result, Eq. (|2.22[) . is only valid for e > e c since a solution for Eq. (|2.8p does not exist for subcritical values 
of the adsorption potential. It is clear, however, that even in the subcritical region, e < e c , there are still monomers 
which occasionally touch the substrate, creating thus single loops at the expense of the tail length. This affects and 
modifies therefore the distribution Ptaii in the vicinity of e c . One can take into account this additional contribution by 
considering a single loop - tail configuration. Pictorially the latter can be inferred from Fig. [4}) where instead of two 
loops and a tail one should imagine a single loop adjacent to the tail. The partition function of such configuration is 



given by Zi-t = 



N — l 
M 3 



Ok 

(jv-i) 1+ * 



On the other side, the partition function of a tail conformation with no loops whatsoever 



(i.e., of a tethered chain) is Z t — /U3 -/V 71 . Thus the probability -P t S[(0 to find a tail of length I next to a single 
loop of length N 



I can be estimated as 



z. 



TV 1-71 



K vw=F* - (2 - 23) 

Evidently, Eq. (|2.23[) predicts a singularity (that is, a steep maximum) in the distribution of tails when I = N. One 
may expect that in the vicinity of the critical point, e w e c , the observed distribution of tails will be given by an 
interpolation between the expressions shown in Eq. (|2.22[) and Eq. (|2.23[) . Hence, the overall tail distribution can be 
represented as 



'^expf-dte-ec) 1 ^ l] 



*W0 = { w 



e > e c 



A 2 N 1 —«i 
l»(N-l) 1 +'t> ' 



(2.24) 



k if> " 



e < e c 



Evidently, close to the CAP this distribution is expected to attain a [/-shaped form with maxima at I w 1 and I = N. 
This shape of Ptail(0 has been predicted earlier for a Gaussian chain by Gorbunov et al.[24|. In close analogy with 
Eq. (|2.24p . the distribution of loops reads 



' _1_ exp [- Cl (e -ec) 1 '* l] 



Ploop(0 = < 



li+<f ' ;i+*(jv-i)' 3 ' 

{ ii+<p(N-iy 



(2.25) 



e < e c 



In Eqs. (|2.24p - (|2.25p A\ , A2 , B\ , B2 are some constants. As we shall see in Section [V] the simulation results for 
Ptaii(0i -Pioop(0 are in good agreement with the predictions, Eqs. (I2.24p - (|2.25p . 
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1. Divergence of the average loop and tail lengths at criticality 

The average loop length is defined by the loop GC-partition function, Eq. (|2.2[) , as 

dlnU(z) 

Li — Z 



Oz 



*lZpllBp (2.26) 



where we have used Eq. (|A2p . Taking into account the polylog function behavior given by Eq. (|Alip with the require- 
ment that 1 < a < 2 as well as the solution for z*, Eq. (|2.12p , one gets 

L**%p°l(±)*S ^oc L_ (2.27) 

where the result Eq. (|2.16p has been used. This result is compatible with the scaling prediction based on Eq. (|2.15p . 
Indeed, close to criticality, L « N/N s . From Eq. (|2.15| ) one obtains then the same result, L oc (e — e c ) 1-1 /*. The free 
energy goes as F = TNhiz* oc —N(e — ec) 1 ^ where one has used Eq. (|2. 12[ ) . On the other hand, the free energy is 
proportional to the number of adsorption blobs, i.e., F oc N/g, where g is the length (number of segments) of the 
blob. The adsorption blobs are defined to contain as many monomers g as necessary to be on the verge of adsorption 
and therefore carry an adsorption energy of the order of ksT each. In result the blob length scales as g oc (e — ec)" 1 /^. 
The size of the adsorbed chain perpendicular to the surface, R±, is nothing but the blob size, that is, R± « g v . Thus 
one obtains 

R± oc K-^r- (2.28) 

(e - e c )"/* 

Consider now the average tail length S. In terms of the GC-partition function for tails, Eq. (|2.4p . it reads 

c d\nQ(z) 

b = z 



dz 



l + $(/3,/i 3 ^) v ' ' 



with the exponent /3 = 1 — 71 = 0.32 < 1. This value of the exponent (3 does not allow for the interaction of the tail 

with other building blocks of the adsorbed chain and will be corrected in Sec lIIBI Using the results, Eqs. (jAlip and 
(I2.12p . the expression for the average tail length Eq. (|2.29p can be recast in the form 

fif«(l-/3)(^)^ 1 oc-J-x (2.30) 

Notably, the exponent (3 drops out of this expression. The corresponding tail size Rs ~ S v scales as 

R * K (73^ (2 ' 31) 

Note that the tail size, Eq. (|2.3ip . scales exactly like the blob (and not the loop!) size, Eq. (I2.28[) . 



B. Role of interacting loops and tails 



As mentioned above, the exponent a, which governs the numbers of loops in the configuration of adsorbed polymer, 
determines also the crossover exponent 4> so that it is of prime importance to know the exact value of a. If the 
surface loops are treated as isolated objects (i.e., loop- loop or loop-tail interactions are ignored), the exponent a — 
1 — 711 = 1.39. Recently Kafri et al. [15[ have shown in the context of DNA melting that the interaction of a loop 
with the rest of the chain increases the loop exponent a. In their work the authors of ref. [l5| essentially used some 
results of the renormalization theory of arbitrary polymer graphs, developed earlier by Duplantier [lgl ]. This approach 
makes it possible to treat also polymer chains which are grafted onto a solid surface. Here we give a short sketch of 
Duplantier's results for a polymer graph located close to the surface and then demonstrate how the loop-loop and 
loop-tail interactions lead to the enhancement of the effective surface loop exponent. 

For an arbitrary self-avoiding polymer graph Q, which is grafted on the surface, it has been shown [l8T | that the 
total number of configurations is given by the standard asymptotic expression: 



(2.32) 
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where N = Mj is the total length of the graph made of AT chains (or edges) of length Mj. The surface exponent 

7 s is given by the following general relationship 

7 s = 1 - u{dC + £ s +V s -l) + Y; ( n L<y k + nlal) (2.33) 

k>l 

where v is the Flory exponent and d stands for the space dimensionality. In Eq. (|2.33|) £ is the total number of 
independent constitutive polymer loops in the graph Q (i.e., the surface loops arc not included in C). C s is the total 
number of extremities of polymer lines upon contact to the surface, and n| are the numbers of bulk and surface 
vertices of order k respectively, thus C s — J2k>i ^ n k- gives the number of surface vertices, i.e., V s = J2k>i n k- 
Finally, and a s k are critical bulk and surface exponents which correspond to the /c-arm vertices. In d < 4 these 
exponents can be calculated analytically via the e-expansion but some of them could be also expressed in terms of 
the conventional exponents v, 7, 71 and 711 (l8| . Figure [3] gives an example of a polymer graph with the specification 
of its topological elements. 

The number of configurations given by Eq. (|2.32j) holds when the lengths of all components M a are large and 
comparable to the total length N. As long as at least one of them becomes small, i.e. M a <C N, then one gets 

Z(S) = ^--. G (^,...,^) (2 , 4) 

where the scaling function G(xi,X2, ■ ■ ■ ,xj\f) has a singularity, provided any of the arguments x a goes to zero. In fact, 
in this limit the polymer graph changes its topology and, therefore, the surface exponent 7 s changes too. In the next 
subsection we show how these results could be used to calculate the effective exponent a which takes into account the 
interaction of a surface loop with the rest of the chain. 




FIG. 3: A polymer graph located close to the surface has the following topological characteristics: C = 1, C 3 = 9, V s = 4, 
ri\ = 3, 713 = 2, nf = 1, n| = 1, n% = 2. One surface vertex is fixed whereas the other vertices may move freely. 



1. Surface loop embedded in an adsorbed chain 

Consider the configurations of a chain (tethered with one end on the surface) in the vicinity of the adsorption 
critical point (see Fig. [4|. Let M be the length of a surface loop while K measures the length of the rest of the chain, 
i.e., M + K = N . The number of configurations of the polymer graph, depicted in FigHp,, is 

Z = »f+ K (M + AT- 1 G (2-35) 

where 7^ is the exponent which could be calculated using Eq.( 12.33]) (see below) and the scaling function G(x) ~ I 
for large M and K. In the case when M/K — > one has a crossover to the polymer graph shown in Fig 0)3 where 
the number of configurations Z ~ ^ (AT) 7 " -1 (1/ 'K') la ~ lh (with 7^ being the surface exponent of the corresponding 
graph). These arguments fix the form of the scaling function which can be written as 

c «*{r : ' :; <**> 

In the case of a small surface loop, embedded in an adsorbed polymer, with M ^> 1 and K ^> 1 (but with M/K <C 1) 
one obtains for the total number of configurations 



M 



7 a -lb 



. - K ll-l ( _ ) „ M M M< - lt M K Rll - K (2 37) 
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FIG. 4: An array of surface loops close to criticality. One of the surface loops of length M in the limit M/N <C 1 is contracted, 
changing the topology of the polymer graph from (a) to (b). This contraction procedure makes it possible to derive the scaling 
function G(x). By similar contraction of a tail the graph goes over from (a) to (c). 



The last result indicates that the total partition function may be factorized to Z ~ Z\ oop Z Test where Z\ oop and Z lesi 
are the partition functions of the small loop and the rest of the chain, respectively. Thus, using the notations of 
Eq. (|2.2p , one obtains f2„ = ^/n < i>~^a ( i.e., the effective exponent a becomes 

oi = 76 -7a- (2-38) 

Now we are in a position to determine the exponents 7* and 7^. Let us assume that the polymer graph in Fig. [4^ 
is made of M subchains (TV — 1 being loops and 1 - a tail). The topological characteristics can be specified as follows: 
C = 0, C s = 2JV — 1, V s = N, ni = 1, n\ = 1, n| = A/" — 1. Earlier it has been shown [l8[ that the critical exponent 
o-f = 2u - 1. With these values Eg. ([233)) yields 

1 s a = 2-M{v + l) + a 1 +a s 1 (2.39) 

The corresponding expression for 7^ can be obtained from Eq. (|2.39[) by the substitution TV — > J\f — 1. This yields 

7b =Z + v-N{v + \) + a 1 +a s 1 (2.40) 

The final expression for the exponent a, given by Eq. (|2.38p . then reads 

« = 7 6 s - 7: = v + 1. (2.41) 

With this theoretical prediction the value of the crossover exponent, given by Eq. (|2.16[) . is determined as: 

(j> = a-l = u= 0.588 (2.42) 

where we have taken the best numerical estimate for the Flory exponent v at d = 3 . A comparison of Eq. (|2.42[) 
with Eq. (|2.17p leads to the important conclusion that, depending on the range of the excluded volume interaction, 
the value of 4> may vary significantly. Indeed, if the interactions affect beads from the same surface loop only then <f) is 
given by Eq. (|2.17[) . otherwise (i.e., when the beads from all loops interact) the value of <j> will be enhanced markedly 
(see Eq. ||232J>). 

One should emphasize, however, that Eq. (|2.41|) does not give an exact value for the exponent a, but rather an upper 
limit only. Indeed, the total number of configurations, given by Eq. (|2.37|) . is estimated by a factorized expression for 
the partition function which takes into account the contribution of a loop and the rest of the chain. As a matter of 
fact this is a Mean Field approach which overestimates interactions at the expense of correlations, reducing thus the 
total number of configurations of a loop. The latter is reflected by an increase of a. The precise value of a therefore 
satisfies the inequality 1 — 711 < a < 1 + v. In the special case of a Gauissian chain both the lower and upper limits 
for a merge while for a phantom chain one has 711 = —0.5 (cf. Section 6.2 in [TH) and v — 0.5. Thus, for Gaussian 
chains one obtains the well known value ef> = 0.5. 

Following the same way of reasoning, one may expect that the exponent for the tail, (3 (see Eq. (|2.29p ). is also 
renormalized due to interaction with the rest of the adsorbed chain. Tail contraction when going from (a) to (c) in 
Fig. [4] enables one to obtain for the renormalized /^-exponent the following relationship 



(2.43) 
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where 7* is the surface exponent of the polymer graph given in Fig. [4J;. Again, if the polymer graph given in Fig. [4^ 
is made of M chains then the exponent 7* for the graph Fig. 2J; becomes 

Yc =i-v-M{v+l) + 2al (2.44) 

Taking into account Eq. (|2.39p . one obtains = 1 — v + <xf — ci, whereby the critical exponents (see [HI) are given 
by 

7-1 

af = v + 7l - 1±1 (2.45) 

The calculation gives finally 

= 71 - 7 + 1 (2.46) 

with 71 « 0.68 and 7 « 1.17 so that (3 « 0.51. As expected, the value of the /3-exponent increases as compared to the 
"isolated tail" case, (3 = 1 — 71 « 0.32. 



2. Comparison with other results 

The result, given by Eq. (|2.42[) , deserves a more detailed discussion. One should point out that, generally, the 
value of <j> for the good solvent case in three dimension has been so far fairly controversial. For example, Monte-Carlo 
(MC) data (albeit for relatively short chains N < 100) on a diamond lattice yield cf> — 0.588 ± 0.03 [25| which is in 
complete agreement with Eq. (|2.42p . A recent MC-investigation (2(| has suggested that the uncertainty in the value of 
<j) might be related to the limited accuracy in the determination of the critical adsorption energy e c . Namely, for the 
bond fluctuation model (BFM), which has been used by Descas, Sommer and Blumen (26)], e c ranges between 0.98 
and 1.01, i.e., within ±2.5%. This relatively small change leads to significant variation of </> between 0.5 and 0.59. 
The same authors have shown that the set of parameters, e c = 1.01 and (j> = 0.59, leads to a more accurate scaling 
prediction. The adsorption of the tethered SAW chain on a simple cubic lattice for chain lengths of up to N = 1000 
(by means of the so-called "scanning method" ) gives: = 0.53 ±0.007 [23. In yet another MC-study, based on the 
pruned-enriched Rosenbluth method (PERM) [28|], it was found that <f> is pretty close to 0.5. However, in a more 
recent study of the same author [2t| one determined for <j> an even smaller value: <j) — 0.484 ±0.002. The value 4> = 0.5 
was also been supported by the MC-simulation results based on the off-lattice model [30]. 

The analytical methods for calculation of <f> are based on the field-theoretical renormalization group (RG) study 
of the semi- infinite n- vector model in the n — > limit. In earlier investigations [3lT l32l . |33| the e-expansion (where 
e = 4 — d) up to order e 2 lead to the prediction (j> = 0.67 which deviates widely from all MC-findings. In a more 
recent investigation the so-called massive field-theory approach at fixed d (i.e., the e-expansion has been avoided) was 
extended to systems with surfaces [34| HH|. The result for the crossover exponent reads (f> « 0.52. Thus we believe 
that the present study elucidates the origin for the diversity of results concerning the precise value of <fi and provides 
a physical background of it. 



III. ADSORPTION UNDER EXTERNAL DETACHING FORCE 



The adsorption of a Gaussian chain on a solid plane under detaching force acting on the chain end has been studied 
first by Skvortsov, Gorbunov and Klushin [36|, [33 in the early 90s. For a Gaussian chain the problem can be solved 
rigorously even for a finite chain length TV. The adsorption-desorption transition is of the first order, however, phase 
coexistence and metastable states are absent. 

Below we apply the GC - ensemble approach to the case of self-avoiding polymer chain adsorption under the presence 
of detaching force. Again, the problem has much in common with the unzipping transition of double-stranded DNA 
[Til ]. When a force / is applied to the free end of the tethered chain, the tail GC partition function in Eq. (|2.ip changes. 
The total GC-partition function is then given by 
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where the tail GC - partition function now takes on the form 

oo 

R(z) = l+J2 Staii(m)z m 

m— 1 

In Eq. (|3.2p Stau(m) is the canonical partition function of the tail under applied force: 

HtauM = % [ d 3 rP m (r) exp(/r ± /T). 
m p J 

Here we take into account that the pulling force is directed perpendicular to the plane (in r^-direction). In Eq. 
P m (r) is the end-to-end distance probability distribution function (PDF) . To estimate this function on large distances 
from the solid plane, i.e., at r±_ ^> R m ss am v (here and in what follows a denotes the length of a Kuhn-segment), we 
assume, following Kreer et al. [38[, that under this condition the PDF is given by the des Cloizeaux expression [4(| 
for the bulk: 



(3.2) 



(3.3) 



r 

R n, 



where the scaling function F(x) is 

F(x) = Bx l exp (-Dx s ) . 
In Eq. (|3.5[) B and D are constants while the exponents 5 and t are given by 

1 



1 - v 



and 



t = 



0-d/2 + dv 
1 - v ' 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



Here = 1 — 71 is the tail surface exponent and d = 3. Note that in the limit r± ;§> R m the only difference between 
the PDFs in the bulk and in the semi-infinite case lies in the fact that instead of the exponent 7 in Eq. (|3.7[) one has 
7i. The integration over the coordinates parallel to the plane in Eq. (|3.3[) is readily carried out and one obtains 



Swi(m) = ^ — / <h-_ [ Tr - 



Rn 



2+t-S 



exp 





00 



-D 



fr± 
T 



dx x 



2+t-S 



exp ( -Dx + f m x 



(3.8) 



where the normalization constant C = <5-D (3+t)/5 ~ 1 /r'[( 3 + _ l \- Tnc integral in Eq. (j3~g| can be tackled 
by the saddle point method (since f m = fR m /T ^> 1). The saddle point itself is defined by the value x sp — 
(fm/ (SD)) 1 /^ -1 ) ~ fm 1 , or, in terms of the r^-variable, 



sp 



Rn 



which is nothing but the well-known Pincus deformation law [4l| . Finally, Eq. (|3.8|) becomes 

exp 



Staii(m) = at (f) e ^ 



rn 



(3.9) 



(3.10) 



with ai and a 2 being constants, the dimensionless force / = fa/T, and the exponent 8 = (2 + 1 — 35/2) /(S — 1). Thus 
the GC-partition function , Eq. ()3.2|) . can be written as 



R(z) = l + ai f e Ap z^eM^f 1 '") 



rn — 1 



l + ai / e $(V,^ 3 exp(a 2 / 1 ^)) 



(3-11) 
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where we have defined the new exponent 

ip = /3-v0=—^--(d-2)u. (3.12) 

One should point out that the exponent drops out from the final expression for ip which for d = 3 is defined as 
if} = 1 — v. 

It is evident from Eq. (|3.11[) that (cf. Eq. (|All[ )) at z — > ii^ 1 cxp(— 02 f 1 ^) the tail GC-partition function has a 
branch point at z — z# , i.e. 



^ r(i-^)(z#)^ 

where 1 — ip < 1 and 



~ oi (/)" \ ~ (3-13) 

z# — z\ 



z* = ^ 1 exp(-a 2 / 1 /-). (3.14) 

Turning back to the total GC-partition function, Eq. (|3.1D . one may conclude that 3(z) has two singularities on the 
real axis Rez: the pole z* which is defined by Eq. (|2.5p . and the branch point z# given by Eq. ()3.14|) . It is well known 
(see, e.g., Sec. 2.4.3 in 20]) that in the thermodynamic limit, TV — > oo, the contribution to the coefficient of z N (i.e., 
to Sjv) consists of contributions by the pole and by the branch singular points, i.e. 

E N ~ d (z*y( N +V + — °* - N~+ {z #)-(N+i-^ (3 . 15) 

r(i-V) 

The singular points, z* and are involved in Eq. (|3.15[) with large negative exponents. Hence, for large N only the 
smallest of these points matters. On the other hand, z* depends on the dimensionless adsorption energy e only (or, 
on w — exp(e)) whereas z# is controlled by the dimensionless external force / (cf., Eq. Q3.14|) ). Therefore, in terms of 
the two control parameters, e and /, the equation 

z*(e) = z*{}) (3.16) 

determines the critical line of transition between the adsorbed phase and the force-induced desorbed phase. In the 
following this line will be referred to as the detachment line. The control! parameters, eo and fu, which satisfy Eq. 
p,16[) . will be named detachment energy and detachment force, respectively. On the detachment line the system 
undergoes a first-order phase transition. The detachment line at fr> — > terminates in the critical adsorption point, 
e c , where the transition becomes of second order. In the vicinity of the critical adsorption point the detachment force 
fo behaves as 

Id ~ (e - e c y^ (3.17) 
where we have used Eq. (|3.16p as well as Eqs. (12.12[) and (|3.14l) . 



A. Order parameter 

Let us study first how the fraction of adsorbed monomers n = N s /N, which we use as an order parameter, depends 
on the pulling force at fixed value of the contact energy t\ > e c . For / < fu it is clear that z* < z# and the first 
term in Eq. p,15[) dominates over the second one. In this case the order parameter 



d hi z*(w) 



dim 



(3.18) 

w— cxp(ei ) 



is constant independent of the force. At / > fr> (i.e., after crossing the detachment line) z* > z# and the second 
term in Eq. (|3.15[) prevails. Since z& is w-independent, it is evident that n = 0, i.e., the polymer is totally detached. 
In result, the n vs. / dependence resembles a step - function with a jump at f = f D , 

Now let us fix the force / = /i and investigate how the order parameter n depends on the adsorption energy e or on 
the fugacity w. Again, Eq. (|3.16p at / = f\ defines a detachment energy e^,. At e < ed one has still z# < z* and the 
second term in Eq. (|3.15[) dominates so that the chain is completely desorbed (i.e., n = 0). At e > erj only the first 
term in Eq. (|3.15p survives so that the relationship n vs. e follows the conventional adsorption dependence without 
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FIG. 5: The 'order parameter' , n, against the: (a) energy of adsorption e (/ is given as a parameter); (b) - against the pulling 
force /, with e a parameter. Vertical lines denote the discontinuous jumps of n, indicating a 1 st — order transition. The noifo) 
(full line) in (b) denotes the order parameter value at the detachment line, (c) The phase diagram of the adsorption-desorption 
transition under puling force / in dimensionless units. An arrow at e c denotes the point of critical adsorption for / = 0. (d) 
The reentrant phase diagram - the same as in (c) but with force against temperature in dimensional units at fixed value of 
e > e c . The largest force / for which chain adsorption may still take place occurs at temperature T max , as indicated by an 
arrow. 



any force-influence. The transition at e = ejj is of first order whereby the order parameter jump grows as the force 
/i increases. 

In Figure^, b we show the predicted variation of the order parameter for an infinitely long chain, following from the 
present consideration. The boundary of the region of adsorbtion, shown in the phase diagram in Figure [SJ;, denotes 
the line of critical values of detachment force for any given attraction of the substrate as described by Eq. (|3. 16|) . 

The adsorption-desorption first order phase transition under pulling force has a clear dichotomic nature (i.e., 
it follows an "either - or" scenario): in the thermodynamic limit N — > oo there is no phase coexistence! The 
configurations are divided into adsorbed and detached (or stretched) dichotomic classes. The metastable states are 
completely absent. Basically, this is in line with the general thermodynamic principles which argue that in thermal 
equilibrium the thermodynamic potentials are convex functions of their order parameters. This exclude multiple 
minima and metastable states [42j. 
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B. Reentrant behavior of the phase diagram 

The results given in Section IIII Al demonstrate that the detachment line on the phase diagram is a monotonous 
function in terms of the dimensionless quantities fu vs. ejj. Recently, it has been revealed that the detachment 
line, when represented in terms of dimensional variables, force fu versus temperature T, goes (at the relatively 
low temperature) through a maximum, that is, the desorption transition shows a reentrant behavior! Below we 
demonstrate that this result follows directly from our theory. 

First, one should note that the low temperature limit implies large values of the ratio e = e/ksT. On the other 
hand, the solution z*(w), which results from Eq. (|2.8p , goes to zero, i.e., z* — > 0, when e — > oo. One may assume 
that under these conditions z*/i2C E — > 1~ (this will be proven a posteriori) . Then the polylog function in the l.h.s. 
of Eq. 1)2.8)1 reads ^(a,fi 3 z*) « fi 3 z* but $ _1 (A, ^wz*) « ci(l - H2WZ*) X ~ X (where we have used Eq. (|A11|) and the 
fact that A < 1). Taking into account Eq. 1)2.8)1 . one arrives at the following result 

H 3 z* kci{1- ^ 2 wz*) 1 - x (3.19) 

This equation determines the function z*{w) at large w. To zero-order approximation the solution reads z? \ ~ 

(/i2w) _1 . Within the first order approximation z*^ ks (jj, 2 w) — S where the decrement 6 is found as S = 

(1/ H2w)(/j,3/ ^w) 1 ^ 1 ^^ . This result is consistent with the assumption z*/i2e e — * 1~ so that the solution of Eq. 
(|3.19[) in the main approximation can be written as 

z* w — e" e (3.20) 

By making use of this solution as well as of the result given by Eq. (|3.14p in the Eq. (|3.16)) . the detachment line at 
large dimensionless detachment energy erj = e/T and force fp = afn/T can be written as 



D 



e D - In 



(3.21) 



Thus, in terms of the dimensionless controll parameters fo increases as the energy en increases. Notably, however, 
if the same detachment line is represented in terms of the dimensional control parameters, detachment force fo vs. 
detachment temperature To (with the dimensional adsorption energy £o being fixed) , one encounters a nonmonotonic 
behavior 



JD — „ > 

a [In \fi 2 

which is shown in Fig. [SJi. The curve given by Eq. (l3.22p goes through a maximum at a temperature given by 



(3.22) 



T max = iL_^l££. ( :-.._):-!! 

In 



(a) 



Such nonmonotonic behavior is termed reentrant and can be observed in the DNA unzipping process [43J, |4J, |45( as 
well as in the case of stretched polymer adsorption on solid surfaces [4^.|4?|. At very low T, however, the expression, 
Eq. (13. 4)1 . for P m (r) [4(| predicts divergent chain deformation [Hj], i.e., it becomes unphysical. One can readily show 
that in this case the correct behavior is given by fa = Sq + T ln(/^3 / /X2) ■ 



C. Average loop and tail lengths close to the detachment line 



As long as the adsorption energy e > e c (or w > w c ), the average loop length L remains finite upon the detachment 
line crossing Namely, at / < fo the fugacity z = z*(w) and the average loop length are given by 



dlnU{z) 



dz 



t *(w) <f>(a,fi 3 z*(w)) 



(3.24) 



Thus, at / < fr> the force does not effects the loop length. At / > fr> the fugacity is given by z — z#(f) where z# 
is determined from Eq. (|3.14|) . In this case the average loop length reads 



dlnU(z) 



dz 



$[a-l,V3Z*{f) 
'(/) $ (a,Ai 3 ^ # (/)) 



(3.25) 
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Recall, that at e > e c and / > /d we have (i^z^ < /13Z* < 1. In this case the function given by Eq. (|3.25|) declines 
when the force grows - see Figure 
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FIG. 6: Variation of the average loop size, L, with detachment force strength / for several values of the adsorption energy e 
(given as a parameter), (b) Mean tail size S against / at different substrate attraction e. 



In contrast, the average tail length S diverges in the vicinity of the detachment line. Indeed, at / < fo the average 
tail length is given by 



S = 



d\xiR{z) 



dz 



z=z*(w) 



aj 9 ® (V - l,M3Z*(u ) )exp(a 2 / 1 ^)) 
1 + ai/ e $ (y,^ 3 z*Hexp(a 2 / lAl ) 



1 -/i 3 z*exp(a 2 / 1 /p) 



(3.26) 



because tp < 1 and fi^z* exp^/ 1 /^) < 1 (cf. Eq. (jAlip ). In the vicinity of the detachment line 1 
^ 3 z*(w) exp(a 2 / 1/A1 ) w {f D - f)/{vf D ) and, therefore, 



a Id 
o oc — -. 

Id- J 



At / > /d the fugacity z — z#(f) and hence, 



S = 



dhxR(z) 



dz 



(3.27) 



(3.28) 



The divergence in Eq. p.28p follows immediately from Eq. (|3.13p which holds in the thermodynamical limit. In practice, 
however, for a large but finite chain length S — » N at / > fr>- Thus, despite the abrupt first order phase transition, 
as far as the order parameter n is concerned, the detachment in terms of the tail length S starts diverging already at 
/ < Id as one comes close to the critical detachment force Jo- 



D. Latent heat variation upon detachment 



What is the internal energy change while crossing the detachment line? At / > Jd the stretching energy E follows 
the Pincus law, so that 



In the adsorbed phase 



E(f = f D +0) = NTpJ" 



E(f = f D -0) = -NT e D n(w D ) 



(3.29) 
(3.30) 
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In result, the latent heat q, consumed upon detachment (or, due to force-induced desorption,) reads 



q = E(f = f D +0)- E(f = f D - 0) = NT 



7l/u 

id 



e D n(w D ) 



> 



(3.31) 



i.e., the heat is absorbed by the system during the force -induced desorption. In the vicinity of the critical point 
f D ~ (e — e c Y^ and n ~ (e — e c ) > thus to a leading order 



g« ATe c (e-ec) 1 / 



(3.32) 



IV. MONTE CARLO SIMULATION MODEL 

We have investigated the force induced desorption of a polymer by means of extensive Monte Carlo simulations. 
We use a coarse grained off-lattice bead-spring model [48] which has proved rather efficient in a number of polymers 
studies so far. The system consists of a single polymer chain tethered at one end to a flat impenetrable structureless 
surface. The effective bonded interaction is described by the FENE (finitely extensible nonlinear elastic) potential: 



Ufene = —K(l - l ) 2 ln 



l-lo 



(4.1) 



with K = 20,l max — l,l Q — 0.7,l min = 0.4. In fact, l max sets the length scale in our model. The nonbonded 
interactions between monomers are described by the Morse potential. 



Um(t) 



exp(-2a(r - r mm )) - 2exp(-a(r - r mm )) 



(4.2) 



with a = 24, r mm = 0.8, e M /k B T = 1. 

The surface interaction is described by a square well potential, 



U w (z) 




(4.3) 



where the range of interaction S = l m ax/^- The strength e of the surface potential is varied from 2.0 to 7.0 and 
k B T = 1. 

We employ periodic boundary conditions in the x — y directions and impenetrable walls in the z direction. The 
lengths of the studied polymer chains are typically 32, 64, and 128. The size of the simulation box was chosen 
appropriately to the chain length, so for example, for a chain length of 128, the box size was 256 x 256 x 256 . All 
simulations were carried out for constant force, that is, in the stress ensemble. A force / was applied to the last 
monomer in the z-direction, i.e., perpendicular to the adsorbing surface. 

The standard Metropolis algorithm was employed to govern the moves with self avoidance automatically incorpo- 
rated in the potentials. In each Monte Carlo update, a monomer was chosen at random and a random displacement 
attempted with Ax, Ay, Az chosen uniformly from the interval —0.5 < Ax, Ay, Az < 0.5. If the last monomer 
was displaced in z direction, there was an energy cost of —fAz due to the pulling force. The transition probability 
for the attempted move was calculated from the change AU of the potential energies before and after the move was 
performed as W = exp(— AU/ksT). As in a standard Metropolis algorithm, the attempted move was accepted, if W 
exceeds a random number uniformly distributed in the interval [0, 1]. 

As a rule, the polymer chains have been originally equilibrated in the MC method for a period of about 5 x 10 5 
MCS after which typically 500 measurement runs were performed, each of length 2 x 10 6 MCS. The equilibration 
period and the length of the run were chosen according to the chain length and the values provided here are for the 
longest chain length. 



V. MONTE CARLO SIMULATION RESULTS 



A. Determination of the detachment point 



In the absence of external pulling force, the transition of a polymer from desorbed to adsorbed state is known to 
be of second order, and the fraction of adsorbed monomers, n, can be identified as an order parameter. Therefore, 



16 




fa/k n T=2.0 



(b) 




Inflexion Point 



Qr].:|Jra::::F:FlCTJ^ 



fa/k R T 



1 2 3 4 5 

e/k B T 



FIG. 7: (a) Plot of the 'order parameter', n, against pulling force /, for an adsorption strength e/ksT = 3.0. The polymer chain 
length is iV=128. The tangent at the inflexion point of the curve meets the abscissa at fn which we define as the detachment 
force, (b) The 'order parameter', n, against the adsorption potential e for fixed pulling force / = 2.0. The tangent at the 
inflexion point of the curve meets the abscissa at Ec/fer where the polymer adsorbs at the surface plane. 



in our computer experiment we use n to determine the point of polymer detachment from the adsorbing surface. At 
constant surface potential, e, one finds that n steeply decreases upon a small increase of the pulling force whereby the 
polymer chain undergoes a transition from an adsorbed phase to a grafted-detached state. In order to locate the point 
of chain detachment, we draw a tangent at the inflexion point of the curve n vs. /. The detachment force, fn, is then 
identified as the point where the tangent intersects the abscissa (/-axis) - see Fig. [2a). Thus one can determine the 
detachment force as a function of the adsorption potential e. Alternately, from the plot of n against the adsorption 
potential e, with the pulling force / held constant, one can observe that as sharp growth of n as the potential is slightly 
increased. The critical potential for chain attachment at the transition point can be found similarly as indicated in 
Fig.J7[b). 

Fig.[8]Ja) shows the variation of the order parameter with changing surface potential for several values of the pulling 
force. Evidently, the larger the pulling force, the stronger the surface potential, needed to keep the polymer adsorbed 
on the plane. In the absence of a force, the order paramenter changes smoothly. For larger forces, however, the 
transition becomes rapidly abrupt. This abrupt behavior of the order parameter is in close agreement with our 
theoretical predictions, depicted in Fig. [5j In Fig. [Sjb) we show the variation of the order parameter n with changing 
force / for various adsorption potentials e. The threshold values for polymer desorption, erj(/) and /13(e), as obtained 
for chains of different length, are then extrapolated to obtain the corresponding values in the thermodynamic limit 
— > 00. 

Our observations show that en increases slightly (i.e., the finite-size effects are rather small) with growing chain 
length N. By extrapolating the data to l/N — > one obtains then en for infinite length of the polymer chain. 
Similarly, the detachment force at fixed surface potential e may be determined in the thermodynamic limit. 



B. Adsorption-desorption phase diagram under pulling 

Using the threshold values of fn and en for critical adsorption/detachment in the thermodynamic limit, one can 
construct the adsorption-desorption phase diagram for a polymer chain. The phase diagram may be obtained by any 
of the two methods, i.e., (i) by fixing of the force and locating en, and/or (ii), by fixing of the surface potential and 
locating the detachment force fn)- The resulting phase diagram is displayed in Figure [5J The inset in Fig [H] shows 
that fn oc (e — e c ) ' 97 which may be compared to the theoretical prediction fn ~ (e — e c ) u ^. Hence, this method 
gives us an estimate for the crossover exponent </>. For e c — 1.67 , we find <fi ~ 0.59 ± 0.02. 
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FIG. 8: (a) The 'order parameter', n, against the surface potential, e, for various pulling forces. The chain has length iV=128. 
(b) Variation of n with the pulling force, /, for several surface potentials. 
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FIG. 9: Plot of the critical detachment force fr> against the surface potential e. In the inset in a double logarithmic plot fo is 
plotted against (e — e c )/fc_gT. The critical adsorption potential for zero force has been found earlier [39| to be e c = 1.67. 

C. Average lengths of loops and tails 

In FigfTUk we plot the PDF of the loop sizes for a chain with N = 128 at several strengths of the adsorption potential 
e in the absence of pulling. One may readily verify that the PDF has a peak for loops of size unity which suggests 
that most frequently single-segment defects (that is, vacancies in the monomer trains) occur in the conformation of 
adsorbed chain. However, for e < e c one may detect clearly in Fig llOa slight increase in the distribution for loops of 
size / N which becomes more pronounced at smaller e ~ 1.0 -j- 1.2 in full agreement with the double-peaked shape, 
predicted by Eq. {222! ■ 

The average loop size L is plotted against the surface potential (with regard to its critical value at the adsorption 
point), (e— e c )/kBT in FigfTUb. We find that, well inside the region of adsorption, L scales as a power law, L oc (e — e) x . 
The exponent x, plotted as a function of N in the inset, is negative, therefore, stronger attraction makes the loops 
smaller while the mean loop size evidently increases with growing chain length N which is a finite size effect. The 
exponent x approaches —0.96 in the limit 1/N — * - see inset in Fig [TUb. This provides another estimate of the 
crossover exponent <j> since x = 1 — l/<fi, according to Eq. (|2.27[) . Thus we find cf> w 0.51 ± 0.02. From Fig 1 10b it is 
evident that the slope of the L vs. (e — e)/kBT curves visibly changes as one comes closer to the CAP. In the immediate 
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FIG. 10: (a) Distribution of loop sizes for chain length N = 128 at different strength of the surface potential and no pulling force. 
The inset shows the same in double logarithmic coordinates. The measured slope at e c /fcsT = 1.67 (full line) is —1.38 ± 0.02 
which practically coincides with the prediction Eq. (|2,18[) . (b) The average loop length plotted against (e — e c )/kBT where 
tc/ksT = 1.67, for various chain lengths in double logarithmic coordinates. The slopes x, indicated by a dashed line, are 
obtained from the L vs. (e — e c )/fcsT curves, and plotted against 1/N in the inset. Extrapolation to 1/N —* yields x ps 0.95. 



vicinity of e c the slope is small and the corresponding estimate for the crossover exponent in this region is <p ~ 0.63. 
One should bear in mind, however, that this is due to the finite length of the chains used in the simulation which 
limits the possibility for the loop size to grow indefinitely, especially at e c . Therefore, we use and depict measurements 
of the slope sufficiently far from the CAP where it tends to a constant value, indicated by the dashed line in FigllOb. 

In Fig. [Tl7a) we plot the PDF of the tail size for a chain with N = 128 at several strengths of the adsorption 
potential in the absence of pulling. An interesting feature of the tail distribution function for e = 1.70 immediately 
at the CAP, e c = 1.67, is the observed bimodal character. It means that there are two dominating chain populations, 
one with few loops and a long tail, and the other with many loops and a very short tail. Our simulation result 
thus confirms the shape of the tail distribution at criticality, Eq. (|2.24p , and appears in excellent agreement with the 
analytic result, derived earlier by Gorbunov et al. [24j |. indicating that in the vicinity of the critical adsorption point 
(CAP) chain conformations are either loop- or tail-dominated. 

In Fig II If b) the average tail length, S, is plotted against (e — e c )/fcsT. Again, S is found to scale as a power law 
with the adhesion strength, S oc (e — e c ) v where y is negative, decreases with N, and approaches eventually —1.67 for 
1/N — > . This result can be compared to Eg. ( 12.30]) . The corresponding estimate of 4> is thus 0.60. 

We turn now to the properties of adsorbed chains in the presence of pulling force. A remarkable feature of the 
probability distribution of the order parameter is the absence of a second peak in the vicinity of the critical strength 
of adsorption, e_o « 6.095 ± 0.03, which still keeps the polymer adsorbed at pulling force fa/ksT = 6.0, see Fig. [T3] 
Somewhat further away from en, one observes a clear maximum in the distribution H(n), indicating a desorbed chain 
with n w 0.01 for e = 6.05, or an almost entirely adsorbed chain with n « 0.99 for e = 6.15. This lack of bimodality 
in the H(n) confirms the dichotomic nature of the desorption transition which rules out phase coexistence. 

In Fig[T3"b, the average loop length, L is plotted against the external pulling force / for e/fc^T = 4.0. For / below 
the detachment threshold, fo, the average loop size appears to be constant independent of the force. As the force 
/ exceeds fz>, the average loop size decreases in close agreement with the theoretical prediction, shown in Fig. [5^,. 
In Fig IT3"b. the average tail length, S, is plotted against the difference fo — f for several chain lengths at surface 
adhesion e/ksT = 4.0 in double logarithmic coordinates. As the applied pulling force / gradually approaches the 
threshold force for detachment, fjj, the tail gets systematically longer and comes close to the length of the chain N. 
Evidently, if one takes into account the finite-size effects which lead to the observed bending of S ks N at stronger 
pulling, the tail S scales as (fo — f)~ w - The exponent w approaches 1.01 (see inset in Fig[T5b) at e/fc^T = 4.0. This 
may be compared to the theoretical prediction of Eq. (|3.27|l which predicts indeed w = 1. 

Eventually, in Fig [T47 a) the PDF of the tail size s is plotted at different strengths of the surface potential e while the 
force, applied to the chain end, is held constant, / = 2.0. In contrast, in Fig [T47 b). we display the distribution of tail 
size s for the case when the adhesion strength is fixed, e/fcgT = 4.0, whereas the pulling force / is varied. Both graphs 
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FIG. 11: (a)Distribution of the tail size for different surface potentials in a polymer of length TV = 128 with no pulling force, (b) 
The average tail length S against (e — e c )/fcsT plotted for various chain lengths in double logarithmic coordinates. The slopes 
obtained from these curves are plotted against 1/TV in the inset and extrapolated to get the thermodynamic limit TV — > oo. 
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FIG. 12: Distribution of the order parameter n for a pulling force fa/ksT = 6.0 and different strength of adhesion ejksT . 
The chain length is TV = 128 and the threshold value of the surface potential for this force is en ~ 6.095 ± 0.03. The values 
e/fcsT = 6.09 and e/ksT = 6.10 are on both sides of the detachment line, cf. Fig. [9] 



are remarkable in that they reflect the transition from fully adsorbed polymer, characterized by a sharp peak in the 
PDF at vanishing tail sizes, to detached chain when the pulling force exceeds the threshold fr> and the corresponding 
PDF is peaked at s/N 1. We emphasize again that although this phase transition of chain detachment is clearly 
of first order, no trace of a bimodal distribution in the vicinity of the transition line can be detected! Thus, the 
states on both sides of the phase boundary /u(e) cannot coexist simultaneously which underlines the peculiar nature 
of this phase transformation. At this point we should like to point out, however, that this exotic feature of the 
detachment transitions has meanwhile been established also in the case of the so called escape transition of a polymer 
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FIG. 13: (a) The average loop length plotted against the pulling force / for fixed e/ksT = 4.0. (b) The average tail length S 
is plotted against the (/ — /d) _1 for various chain lengths for e/ksT = 4.0 in double logarithmic coordinates. The inset shows 
the extrapolated slope for N —* oo go to unity, as predicted by Eq. (|3.27p 




coil, deformed under the tip of an Atomic Force Microscope [49|, It has been shown rigirously recently [5l|, that 
despite its first order nature, the escape transition takes place without phase coexistence. Most probably, this unusual 
feature is due to the topological connectivity of polymer chain as quasi-onedimensional systems. 



VI. SUMMARY AND DISCUSSION 



In the present investigation we have studied the force-induced desoption transition of a polymer chain in contact 
with an adhesive surface. We treat the problem within the framework of the Grand Canonical Ensemble approach and 
derive analytic expressions for the various conformational building blocks, characterizing the structure of an adsorbed 
linear polymer chain, subject to pulling force of fixed strength. Closed analytic expressions for the fraction of adsorbed 
segments (i.e., the order parameter of the desorption transition) and for probability distributions of trains, loops and 
tails have been derived along with expressions for the corresponding first moments in terms of the surface potential 
intensity both with and without external force. As expected, all these conformational properties and their variation 
with the proximity to the CAP are governed by a crossover exponent <f>. 
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A central result in the present work is the calculation of cj) using the approach of Kafri et al. [l5| which provides 
insight into the background of the existing controversial reports about its numeric value. We demonstrate that the 
value of cj> may vary within the interval 0.39 < 4> < 0.6, depending on the possibility of a single loop to interact with 
the neighboring loops in the adsorbed polymer. Since this range is model-dependent, one should not be surprised 
that different models produce different estimates of <f> in this interval. 

A comparison with the results from extensive Monte Carlo simulations demonstrates the good agreement between 
theoretic predictions and simulation data. 

In particular, we verify the gradual transition of the PDF of loops from power-law to exponential decay as one 
moves away from the critical adsorption point to stronger adsorption. We demonstrate that for vanishing pulling 
force, / — * 0, the mean loop size, L cx (e — e c ) 1- *, and the mean tail size, S cx (e — e c )*j diverge when one comes 
close to the CAP. In contrast, for a non-zero pulling force, / ^ 0, we show that the loops on the average get smaller 
with growing force while close to the detachment threshold, / w fr>, the tail length diverges asS'oc(l — 

Eventually, we derive the overall phase diagram of the force-induced desorption transition for a linear self-avoiding 
polymer chain and demonstrate its reentrant character when plotted in terms of detachment force fr> against system 
temperature T. We find that despite being of first order, the force-induced phase transition of polymer desorption 
is dichotomic in its nature, that is, no phase coexistence and no metastable states exist. This unusual feature of 
the phase transformation is unambiguously supported by our simulation data, e.g., through the comparison of the 
order parameter probability distributions on both sides in the immediate vicinity of the detachment line whereby no 
double-peaked structure is detected. 

Finally, we should like to to emphasize that while the present investigation will hopefully shed new light on the 
force-induced desorption transition of a linear polymer from a sticky surface, a lot more work is needed before a 
comprehensive understanding of this phenomenon is achieved. In this work simulations have been carried out within 
the framework of a constant force ensemble. In their comprehensive treatment of the problem, however, Skvortsov et al. 
[HI have shown that one may well work in the constant height ensemble whereby one uses the end-monomer /i-position 
as an independent parameter and measures the force, exerted by the chain on the end monomer. Notwithstanding 
the equivalence of both ensembles, some quantities behave differently in each ensemble and this becomes evident 
only if one presents the salient features of the system behavior in each particular ensemble. Thus in the fixed-height 
ensemble, a different and rather interesting thermodynamic behavior of the measured mean detachment force and 
of the fraction of adsorbed segments agains h is expected to be observed. Typically one then observes a constant 
force plateau while the height of the chain end monomer is varied. Such a behavior can be inferred even within the 
fixed-force ensemble. Indeed, as shown in Fig. [Mb the PDF P(s) is practically flat for the critical detachment force 
Id = 3.40, meaning that all chain end heights are equally probable at this particular force. The latter is equivalent 
to a constant-force plateau in the fixed-height ensemble. A verification by computer experiment is among our tasks 
in the immediate future as well as a study of the so far unexplored kinetics of chain detachment. 
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APPENDIX A: PROPERTIES OF THE POLYLOG FUNCTION 



The poly log function $(a, z) is defined by the series 




(Al) 



k=l 



which converges at \z\ < 1. From the definition, Eq. (|A1|) . one immediately obtains 

z <&(a, z) = $(a - 1, z) 
az 

The calculation of the series Eq. (|A1|) (see Sec. 1.11 in ref. [2lJ) gives 



(A2) 




(A3) 
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where T(x) is the gamma- function, £(x) is the Riemann zeta-function, and the exponent a is noninteger, i.e. a ^ 
1,2,3,... 

Consider now the case of integer values of a. The gumma-function r(x) has poles at all negative integer arguments 
whereas the pole of ((x) 1S placed at x — 1. One may write a — to + 6 where m is a positive integer and 5 — ► 0. Then 
in the vicinity of the poles the gamma- and zeta-functions can be rewritten as 



T(l - m - 5) = 



((1 + 5) = 



(~l) ro 
(to - 1)! 



(A4) 



where ip(x) is the digamma function (or ^-function) defined as the logarithmic derivative of the gamma-function, 
ip(x) — dlnT(x)/dx. One should also take into account that 
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= 1 + Shx 











0(6) 



(A5) 



After taking into account Eqs.jA4| and (|A5|1 in Eq. (|A3|) and due to the cancellation of poles in the gamma- and 
zeta-functions at small values of S the polylog function, Eq. (|Alj) becomes [2lJ 
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r=0 



where the prime indicates that the term r = to — 1 is to be omitted. 

We are interested in the behavior of <&(«, z) at z — > 1. In this case ln(l/z) = — ln[l — (1 — z)\ w (1 — z). At a < 1, 
the main contribution comes from the first term in Eq. (IA3|) . i.e. 
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r(i - a) 
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At a = 1 and z — > 1, and making use of Eq. (|A6[) . one obtains 
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(A7) 
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Finally, at a > 1 the polylog function ^(q!, z) has no singularity at z 
expansion 



1 and Eq. (|A3[) results in the following 



$(a, z) « C(«) + r (! - «)(! - Z T^ - C(« - 1)(1 - 2) + 



(A9) 



In a bit more specific case when 1 < a < 2 we will use the well known relationship T(l — a) = —ir/\T(a)\ sin(7ra)|] so 
that 



$(a,z) « ((a) 



-(i-zY 



C(a -!)(!-«) + 



r(a)| sin(7ra)| 

Taking into account the Eqs. (|A7j) . (IA8|) and (|A10|) . the expression for the polylog function at z — > 1 reads 

at a < 1 



(A10) 



f r(l-a) 

(1- 2 )1— ' 



X(a) - a a (l - z)*- 1 - b a (l - z) + 
where the coefficients a a = 7r/r(a)| sin(7ro!)| and b a = ((a — 1). 



at a = 1 

at 1 < a < 2 



(All) 
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